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ABSTRACT 



We have studied the effect of gravitational lensing on the Cosmic Microwave 
Background (CMB) anisotropy in flat and open universes. We develop a 
formalism to calculate the changes on the radiation power spectrum induced by 
lensing in the Newtonian and synchronous-comoving gauges. The previously 
considered negligible contribution to the CMB radiation power spectrum of the 
anisotropic term of the lensing correlation is shown to be appreciable. However, 
considering the nonlinear evolution of the matter power spectrum produces 
only slight differences on the results based on linear evolution. The general 
conclusion for flat as well as open universes is that lensing slightly smoothes 
the radiation power spectrum. For a given range of multipoles the effect of 
lensing increases with Q but for the same acoustic peak it decreases with Q. The 
maximum contribution of lensing to the radiation power spectrum for / < 2000 
is ~ 5% for Q values in the range 0.1 — 1. 

Subject headings: cosmic microwave background - gravitational lensing - 
large-scale structure of universe 
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1. Introduction 

Cosmic Microwave Background (CMB) temperature anisotropies, detected for the 
first time by Smoot et al. (1992) with the COBE-DMR experiment, are beheved to be 
generated by the interaction of matter density perturbations and radiation to first order 
in perturbation theory. Numerical codes used to solve the linearized Einstein-Boltzmann 
coupled equations are able to calculate the radiation power spectrum with an accuracy 
better than 1% (See e.g. Sugiyama 1996, Seljak & Zaldarriaga 1996, Bond 1995). Nonlinear 
density perturbations make a small contribution through the Rees-Sciama effect which, 
except for the case of reionization, can be constrained to be <^ 1% (Martinez-Gonzalez, Sanz 
and Silk 1992, Sanz et al. 1996, Seljak 1996a, Tuluie, Laguna & Anninos 1996). However, 
the effect of gravitational lensing on the CMB anisotropies, not included in the numerical 
codes, may appreciably affect the radiation power spectrum. 

Many groups have studied the lensing of the microwave photons using different 
analytical and numerical approaches (Blanchard & Schneider 1987; Cole & Efstathiou 1989; 
Sasaki 1989; Tomita & Watanabe 1989; hinder 1990a, b; Cayon, Martinez-Gonzalez & Sanz 
1993a, b; Fukushige, Makino & Ebisuzaki 1994; Seljak 1996b). They arrive at different 
conclusions about the importance of the effect: the result depends on the particular 
cosmological model considered and on the asumptions made in the calculation. Cayon et 
al. 1993a,b present the formalism to obtain the lensing of the microwave photons by the 
large scale matter distribution in a flat universe with null/non-null cosmological constant. 
However they erroneously used the photon deflection angle instead of the photon angular 
excursion on the last scattering surface relative to its observed value, which leads to a factor 
of a few overestimate of the relative dispersion between two photons (Seljak 1996b, Muhoz 
& Portilla 1996). Some of the previous studies have used models that may not be a realistic 
representation of the large-scale structure observed (e.g. the models used in Fukushige et 
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al. 1994). Another relevant ingredient of the calculations is to appropriately account for the 
evolution of matter density perturbations. Recently, Seljak (1996b) has done a relevant step 
in solving those shortcomings of previous studies. Based on a power spectrum approach 
he includes linear and nonlinear regimes of the matter evolution in realistic cosmological 
models and generalizes the formalism to open universes. However, results on the radiation 
power spectrum are not presented for open universes. Moreover, the nonlinear power 
spectrum evolution considered in that paper is not valid for spectral indexes n < — 1 (as in 
the case of CDM for small scales) and for f] < 1 universes (Peacock and Dodds 1996). 

In this paper we present a formalism to calculate the lensing effect in flat and open 
cosmological models and in two different gauges. Except for velocity and acceleration terms 
associated to the observer and the source which either do not contribute or the contribution 
is negligible, we show that the equations which provide the lensing effect are the same 
for the conformal Newtonian and synchronous-comoving gauges. Results for the effect of 
lensing on the radiation power spectrum are presented for CDM models with 0.1 < ^2 < 1. 
We consider linear and nonlinear evolution for the matter power spectrum. The structure of 
the paper is as follows: in section II we describe the formalism to calculate the gravitational 
lensing effect. The results obtained for CDM open models are presented in section III. 
Finally, the main conclusions are related in section IV. 

2. Formalism 

2.1. Geodesies in the conformal Newtonian gauge 

We will consider the propagation of photons from recombination to the present time, 
the universe being a perturbed Friedmann model with a dust {p = 0) matter content. We 
shall not consider a cosmological A-term, but the generalization to include A 7^ is very 



-5- 



easy. For scalar perturbations, the metric in the conformal Newtonian gauge is given in 
terms of a single potential 0(r, x) as follows 

we take units such that c = SnG = ao = 211'^ = 1 and k/{A \ 1 — ^2 |) = 0, —1, +1 denote 
the flat, open and closed Friedmann background universe. The gravitational potential 
satisfies the Poisson equation 

( V + 3k)<j) = ^p,a^S, (2) 

where 6 is the density perturbation. The Green's function associated to the previous 
equation can be found in the literature (D'Eath 1976, Traschen and Eardley 1986). We are 
interested in the effect of gravitational lensing on high multipoles (/ ~ 10^) of the CMB. 
Only the smaller scales are contributing to such effect, so curvature will show related to the 
angular distance. In fact, the Green's function on such scales can be approximated by 

G{x,x) c^- — \x-x\^ , (3) 

where the distance between the two points x = An, x = X n (being n and n two unit 
vectors in the directions of observation) is given by the equation (A ~ A ) 

\x — X \c — Is"^ + s"^ — 2ss cosaY'"^, s = ; ^, ^„ , cosa = n-n. (4) 

^' 1 — (1 — iZ)A^ 

On the other hand, after a straightforward calculation, the geodesic equation associated to 
the metric (1) gives the following equation for the vector s^ = ^ = ^ 



ds' 
d7 



1 



2 



dt 



k-f-^ {x ■ s)s' - --f^x' - 2k-f(f)x' + 2 -- + 2{V(f)-s) s'- 2-f\V(l)y. (5) 



dr 



Assuming a perturbation scheme ("weak lensing"), this equation can be integrated in the 
form 

X = \n + e, (6) 



-6- 



where n is the direction of observation and A is the distance to the photon for the 
background metric, i.e. 

A = ro-r (A; = 0), X = {I - nyHanh[{l - n){To - t)] {k = -I). (7) 

The perturbation e can be decomposed in a term parallel to n and a term orthogonal 
to such a direction d±. The last term satisfies the following differential equation when 
parametrized by A 



cPa± k [ ^da± 



+ 



-A—— + a_L 
a A 



-2V 



«(A) = ;, ,j>. ^ w{\,\') = {\-\') j_''v. - (10) 



rfA2 27 

where (V_l0)* = (t^*-' — n^n^)d(f)/dnK The solution to the previous equation with the initial 
conditions: a^(A = 0) = = ^(A = 0) is 

5^ = -2 / dX'W{X, A')V±0(A', f = X'n) (9) 

Jo 

where W{X, A ) is a window function 

^^#^, 1^(A,A')HA-A')i±^, 

For photons that are propagated from recombination, A^ = [1 — (1 + r22;r)""'^/^][l — (1 — 
r2)(l + VlZj.)^^^'^]^^ with Zr ~ 10^, to the observer, Aq = 0, the lensing vector /3 is defined in 
the usual way (see Figure 1) 

P = n- 5~5| , (11) 

SO we find 13 = —j-a^{Xr) and the final result, taking into account equations (9, 10), is 

^=-2£dXW{X)V^<P{X,x = Xn) W{X) = {1 - X)^^^^^^, (12) 

because A^ — 1 for Qzr > 10^. 

The lensing vector for a flat universe has been given by Kaiser (1992). For the open 
case, Pyne & Birkinshaw (1996) and Seljak (1996b) have used a window function W that 
agrees after a straightforward calculation with our equation (12). 
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2.2. Geodesies in the conformal synchronous-comoving gauge 

Once we have obtained the expression for the trajectory of the photon in the conformal 
Newtonian gauge, it is easy to calculate everything in the conformal synchronous-comoving 
gauge. The infinitesimal transformation connecting both gauges is 

r' = r + e%T,x), x'' = x' + e\T,x), 6° = ^^(a0), e* = ^V(a0). (13) 

a-^Pb OT a-^pb 

The expressions for e" and e* can be obtained taking into account that in the synchronous- 
comoving gauge one has zero velocity (f* = 0), i. e. 

«o = (1 - e°)M°', u' = -e'u"', v' = -e\ ' = ^, (14) 

and the metric has the following components: qq'o' = — a^(r'), (^o'i' = 0. So, 



a 
— ( 
a 



'' ' V = 0, e°.-7-% = (15) 



and integrating the last equations we get the result mentioned above for (e°, e*). Moreover, 
the metric in the conformal synchronous-comoving gauge reads 



ds^ = a\T'){-dT'^ + 7-2 [(1 - 20 - 2-e° + k-f^^x ■ e)6ij - e^j - ej^^]dx'' dx^'}. (16) 



This last expression for /c = agrees with the one given by Sachs & Wolfe (1967). By 
changing the gauge, the new lensing vector f3 is given by an equation similar to (11), so we 
obtain 

^- = /5-i(e-x..-?x,)-(fX. (17) 

where e*^, e*are given by equation (13). The velocity of the fluid in the conformal Newtonian 
gauge is given by -u = — i^, from which e*= —{a/df)v (/ = dlnD/d\na, D{a) being the 
growing mode). Taking this into account one can easily understand that the new terms 
appearing in equation (17) can be interpreted as Doppler contributions at recombination 



and at the observer and an acceleration term at the observer. For a flat model (A; = 0), we 
explicitly have 



/?-/5-i 



Vr - Vn 



where the linear gravitational potential, 0(x), is time-independent and Vo = — |V0o, 
Vr = — 1(1 + -2r)^"'^V0r and Cq = (^)o- The ratio of these terms, as they appear in equation 
(18), to the angular scale is negligible {vo is given by the Doppler velocity respect to the 
CMB and So can be estimated from our local infall towards either the Virgo cluster or the 
Great Attractor). A similar reasoning can be applied to open universes. Therefore, the 
lensing vector (3 in the synchronous-comoving gauge (that is the appropriate one from the 
point of view of the observations) is approximately given by j3, as defined by equation (12). 

2.3. The influence of weak gravitational lensing on the Ci s 

The correlation function C{6) including gravitational lensing is calculated as the 
average 

€{9) = (A(n + /3(n)) A(n' + /3(n ))), (19) 

where A(n) is the temperature anisotropy field, n and n are two directions such that 
n ■ n = cos 6. By introducing 2D-Fourier components of the temperature anisotropics A^ 
and assuming that the anisotropics. A, and the lensing vector, j3 are uncorrelated, we obtain 

0(9) = ^j^ dqqP^{q){Uqv)), (20) 

where Jq is the Bessel function, z/ = \n — n + j3{n) — (3{n )| and Pa is the 2D-power 
spectrum of the radiation field: {A^AZi) = -Pa(q')^^(^^ Q )• On the other hand, assuming 
weak gravitational lensing, i.e. on the average the relative lensing vector is very small as 
compared to the angle 6, we can make a series expansion in the previous equation obtaining 

C{e) - 0(9) - ^|[Q. - Q^^-gr\^^f- + [Q'^^-0^\^ —^]- (21) 
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Qij is the bending correlation matrix 

Q,, ^ ([A(n) - A(n )][/3,(n) - /3,(n')]), 
and can be decomposed into the trace and an anisotropic component 



(22) 



Ql^2a\e), m^Q 



'' e^ 



— a 



(23) 



where a {9) is the bending dispersion and ^{9) is the anisotropic correlation {^{9) corresponds 
to Cgi^2{9) in Seljak 1996b). Therefore equation (21) can be rewritten as 



C{9) - C{9) 



a 



2 .^2 



d^C{9) ldC{9)i i\(fC{9) ldC{9) 



d9^ 9 d9 



+ 



d9^ 9 d9 



(24) 



On the other hand, taking into account the expansion 



C{9) = —Y,^{2l + l)CiPi{cosi 



(25) 



and the approximation Pi{cos9) ~ Jq{19) for / ^ 1, we get 

Ci^2ti I d99C{9)jQ{l9) (/ > 1). 
Jo 

From equations (24) and (26) 
Ci-Ci = -^Er(2^' + 1)^''^/' [ d99.Ml9)W\0).Ml'9) - mUl'0)i (/ » 1). 



(26) 



(27) 



The next step is the calculation of the dispersion and correlation of the lensing vector 
as a function of the power spectrum P{a, k) defined by 



{6,:{a)6l,{a)) = P{a,k)6%k - k'). 



(28) 



From equation (12) one can obtain 



{(3,(3,) = 4D,D, f rfA^ / rfA'^^Q(A,A',r), 
JO A JO A 



(29) 
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where Di = {Sf — n^ni)d / dn^ and C^{\,\ ,r) is the correlation of the gravitational 
potential at two different times. If one assumes Limber's approximation (see also Kaiser 
1992), i.e. only a small region r with A ~ A is contributing, the previous equation can be 
approximated by 

'•1 rW^(A)i2 roo 



{(3i(3,) = SDiDj / d\ 



A 



es 



drrij-' - e's'y''^[i - (1 - n)y]c^{\ r). (so) 



Notice that the correlation depends only on a single time and s is given by equation (4). 
Introducing the power spectrum P(f,{a, k), the last expresion becomes 



{(3,(3,) = -D,D, f' d\ 

71 Jo 



W{X) 



A 



[l-{l-n)\^] dkk-^PJ\,k)Jo{ks9). (31) 
Jo 



P^{a, k) is given by the Poisson equation (2) for scales (A;^ ^ 12(1 — Vt)) 

P^{a,k)c^{—) k-^P{a,k), 



(32) 



where P{a, k) is the power spectrum associated to the matter perturbations. In the linear 
regime: P{a, k) = D[a)P{k), D[a) being the growing mode normalized to the present time 
(see Peebles 1980). 

Finally, calculating the derivatives that appear in equation (31) and applying equations 
(22, 23) one can obtain 



,^(„=i?^r*rdA 



71 Jo k Jo 



W{\) 



P{a,k) 



l-(l-fi)A2 



l-Jo + isin^^Jo-sin^-Js], (33) 



e(^)^.^+^r^/^Ar^(^) 



71 Jo k Jo 



P{a,k) 



l-{l-n)X^ 



[cos ^ - 3)(1 - Jo - J2) - sin^^Jo 



(34) 



where the argument of the Bessel functions Jq and J2 is ks6. Notice that the behaviour of 
a (9) and ^^^'^{9) for small 9 is linear: 



a{9) 



a9 



e/\9) 



b9 , 6 ~ ^ , 
V2 



9 ' 9 

as will be shown by the numerical calculations presented in the next section. 



(35) 
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3. Results 



With the formahsm presented in the previous section, we have calculated the dispersion 
of lensing a(9) and the anisotropic term of the correlation of lensing ^^^"^(9) as given 
by equations (33-34). We assume a CDM model with a primordial Harrison-Zeldovich 
spectrum, a Hubble parameter h = 0.5 {H = 100/ikm s~^ Mpc~^) and fiat and open 
universe models. The radiation power spectrum not including lensing is normalized to the 
2-year COBE-DMR map as given by the analysis of Cayon et al. (1996) (this normalization 
does not appreciably change with the 4- year data). However, since the lensing effect is 
generated by small scales, << 100 Mpc, it might be more sensible to use the normalization 
(Tg = 0.6, 1, 1.4 for universes with Q = 1,0.3,0.1 following Viana and Liddle (1996). This 
normalization is based on the cluster abundance (see also White, Efstathiou and Frenk 
1993, Eke, Cole and Frenk 1996). For the nonlinear evolution of the power spectrum we use 
the recently improved fitting formula given by Peacock and Dodds (1996). That formula is 
based on the Hamilton et al. (1991) scaling procedure to describe the transition between 
linear and nonlinear regimes. It accounts for the correction introduced by Jain, Mo and 
White (1995) for spectra with n <^ —1 and applies to fiat as well as to open universes. 

In figure 2 it is shown the relative dispersion a (9)/ 9 and the anisotropic term ^^^'^{9)/9 
for three values of the density parameter Q = 1,0.3,0.1. Linear and nonlinear matter 
evolutions have been considered for comparison. Discrepances between the two regimes can 
be noticed at scales 9 ^ 3'. At scales 9 ;^ 6" ct{9)/9 as well as (^'"^{9)19 are below 20% 
being slightly larger as Vt increases. Also, notice that {^^^"^{9)19)0^0 ~ -j^{a{9)/ 9)0^0 in 
all cases (as expected from the considerations made in the previous section). Therefore, 
the anisotropic term should in principle be considered when calculating the distortions on 
the radiation power spectrum contrary to the isotropic approximation often made in the 
literature (we confirm this statement below). Lensing becomes negligible at angular scales 
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above a few degrees. 

The radiation power spectrum including and not including lensing is given in figure 
3a. It is clear from this picture that the effect of lensing is to slightly smooth the main 
features appearing in the spectrum, in particular the secondary acoustic peaks (also called 
Dopplcr or Sakharov). The relative changes of the spectrum due to lensing as a function 
of the multipole / are shown in figures 3b, c. The changes produced grow with / and in 
the same range of / increase with Q. Note, however, that for the same acoustic peak the 
variation increases when Q decreases and the reason for this is the smaller scales involved 
for which the lensing effect is more effective. Considering nonlinear evolution does not 
change the lensing contribution to the Ci for / <^ 2000 as can be seen in figure 3b. We have 
also computed the contribution of the isotropic term to the Ci coefficients and the result is 
shown in figure 3c. This contribution is slightly smaller than the total effect and for some 
multipoles the discrepancy can be significant. Therefore, in general both terms, isotropic 
and anisotropic, should be considered in the calculation of the radiation power spectrum 
with lensing. 

The effect of lensing can be as much as ~ 2% for multipoles / <^ 1000 and ^ 5% 
for / ^ 2000. Therefore, if one wants to compute the radiation power spectrum for a 
particular cosmological model with an accuracy better than 1% such effect should be 
considered. Bending of the microwave photons due to the large-scale structure should be 
considered when analysing data provided by future very sensitive CMB experiments (e.g. 
COBRAS/SAMBA). 
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4. Conclusions 

A formalism has been developed to calculate the lensing effect on the primary CMB 
radiation power spectrum. This formalism provides an expression for the lensing vector in 
fiat and open universes which is approximately the same in both the conformal Newtonian 
and comoving-synchronous gauges. In particular, we give the window function W for open 
models in terms of the distance to the photon from the observer. 

The influence of gravitational lensing on the Ci s has been obtained in terms of the 
bending dispersion cr and anisotropic bending correlation ^. It is found that the contribution 
of ^ to the lensing distortion of the radiation power spectrum is smaller than that of a. 
However this contribution is not negligible and should be considered in the calculation of 
the radiation power spectrum with lensing. 

We use the recently improved fitting formula for the evolution of the nonlinear matter 
power spectrum which provides an accuracy better than 12% for the scales considered 
(Peacock and Dodds 1996). This improvement over previous works (Peacock and Dodds 
1993; Jain, Mo and White 1995) generates a larger lensing dispersion at small scales for 
open models, as compared with Seljak (1996). In spite of this, the contribution of nonlinear 
evolution to the distortion of the radiation power spectrum is negligible. 

For fiat as well as open universes, the effect of lensing is to slightly smooth the primary 
radiation power spectrum of the CMB. For a given range of multipoles the relative change 
of Ci due to lensing increases with Q. However, for the same acoustic peak it decreases with 
Q. The maximum contribution of lensing to the radiation power spectrum for / < 2000 is 
~ 5% for Q values in the range 0.1 — 1. Therefore, the effect of lensing should be considered 
in analyses of CMB anisotropy data provided by future very sensitive experiments. 

We would like to thank N. Sugiyama for providing us with the radiation power 
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Fig. 1. — Diagram describing the geometry and the angles involved in lensing calculations 

on the CMB. Note that e = {xr — Xo)/\xr — Xo\ as appears in equation (11). 

Fig. 2. — Ratios of the bending dispersion (a) and the anisotropic correlation (b) to the 
angular distance, a{9)/9, ^^^'^{9)/9, as a function of the angular distance 9 in units of 
arcmin. Solid, dashed and dotted lines correspond to i7 = 1,0.3,0.1 respectively. Thick 
curves represent the results considering nonlinear evolution whereas thin ones outline the 
results from linear evolution. 

Fig. 3. — (a) Radiation power spectrum including (solid) and not including (dashed) lensing 
for Q = 1,0.3,0.1. For / = 1000 the lower curve represents to 1) = 1 and the upper one 
to fi = 0.1. Relative change in the radiation power spectrum due to lensing for f] = 1 
(solid), Q = 0.3 (dashed) and Q = 0.1 (dotted), (b) Thick curves represent the total effect 
(including nonlinear evolution and the anisotropic term) whereas thin ones outline the result 
of not considering nonlinear evolution, (c) Thick curves represent the total effect whereas 
the thin ones outline the result of not considering the anisotropic term. 
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